Resonance approximation and charge loading/unloading in adiabatic quantum 
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Quantum pumping through mesoscopic quantum dots is known to be enhanced by resonant trans- 
mission. The pumped charge is close to an integer number of electrons when the pumping contour 
surrounds a resonance, but the transmission remains small on the contour. For non-interacting elec- 
trons, we give a quantitative account of the detailed exchange of electrons between the dot and the 
leads (to the electron reservoirs) during a pumping cycle. Near isolated distinct resonances, we use 
approximate Breit-Wigner expressions for the dot's Green function to discuss the loading/unloading 
picture of the pumping: the fractional charge exchanged between the dot and each lead through a 
single resonance point is related to the relative couplings of the dot and the leads at this resonance. 
If each resonance point along the pumping contour is dominated by the coupling to a single lead 
(which also implies a very small transmission), then the crossing of each such resonance results in 
a single electron exchange between the dot and that lead, ending up with a net quantized charge. 
When the resonance approximation is valid, the fractional charges can also be extracted from the 
peaks of the transmissions between the various leads. 

PACS numbers: 73.23.-b, 73.63.Rt, 72.10.-d, 73.40.Ei 
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I. INTRODUCTION 

There has been much recent experimental-' 2 ' 3 and 
theoretical 4 . 5 ! 6 . 7 . 8 . 9 . 10 ! 11 . 12 ! 13 ! 14 ! 15 ! 16 ! 17 interest in adia- 
batic quantum pumping through mesoscopic electronic 
devices, such as quantum channels or quantum dots 
(QDs). Typically, the QD is connected via leads to 
several electron reservoirs, and is subject to a slowly 
varying oscillating potential, with period T = 2ir/u>. 
Under appropriate conditions, the device yields a non- 
zero dc time-averaged current between pairs of termi- 
nals, even when the terminals have the same chemical 
potential. Under ideal conditions, the charge Q trans- 
ferred between the terminals during a period T may be 
'quantized', i.e. very close to an integer times the elec- 
tron charge e. Several recent theoretical studies have 
considered enhancement of the adiabatic pumping cur- 
rent due to resonant transmission through the QD, both 
for non-interactinga*ii*iSii2*iS and interacting electrons 2 *. 
Connections between pumped charge quantization and 
resonant transmission have been reported in different 
contexta 11 . 16 . 18 . 19 . 21 . 22 . 

Usually, the oscillating potential is characterized by 
several time dependent parameters, {Xi(t)}. As time 
evolves during one period T, these parameters follow a 
closed contour in the parameter space. A schematic ex- 
ample is shown in Fig. ^ for two such parameters. In 
parallel to discussing pumping, one can also consider the 
conductance between pairs of terminals, generated by an 
appropriate bias. This conductance, which depends on 
the parameters {A^}, may have resonance peaks in the 
same parameter space. In this context, one freezes the 
time dependence, and considers the conductance at some 
instantaneous values of the {Ai}'s. It has been arguedii 
that the pumped charge Q will be close to being quan- 
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FIG. 1: (color online) Schematic picture of a two-dimensional 
pumping contour, crossing the resonance line at two resonance 
points (B and D). The transmission is maximal at the point 
M. 



tized if the pumping contour surrounds such a peak (e.g. 
at the point M in Fig. ^) , while staying at points with a 
low conductance. 

In the present paper we present an approximate theory 
for adiabatic pumping of coherent non-interacting spin- 
less electrons, which is valid for discrete and distinct res- 
onances, and use this approximation to obtain physical 
insight into the reasons for this quantization. Given a 
conductance peak (e.g. at the point M in Fig. QJ, one 
can usually also identify a 'resonance line', along which 
the conductance decreases from its peak more slowly than 
along other directionsii* 1 ^. Such a line is illustrated by 
the dashed line in Fig. ^ In the example shown in this 
figure, the resonance line is crossed by the pumping con- 
tour twice, at points B and D. Measuring the instanta- 
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neous biased conductance between the two relevant ter- 
minals for each time t during the oscillation period, one 
expects two local peaks at these two 'resonance points'. 
Under appropriate conditions, which include the limit of 
weak QD-terminal coupling, most of the pumped current 
arises when the parameters are close to these resonance 
points: for example, one can identify a 'loading' of the 
QD by some charge AQ r ° s , coming from terminal a, at 
the point B, and an 'unloading' of the QD, by AQ™, S , 
into terminal a', at the point D. The resulting total 
pumped charge per period approaches a robust, detail- 
independent value Q R , which is determined only by the 
ratios of the coupling strengths between the QD and the 
different reservoirs at the resonance points. We also show 
that Q R can be related quantitatively to the measured 
values of the peak conductances. Q R is (almost) quan- 
tized (in units of e) when there is one dominant coupling 
for each resonance. 

Our results can be summarized in a very simple and 
physically transparent way, by considering the occupa- 
tion numbers of the quasi-bound state on the QD, corre- 
sponding to each transmission resonance. Each time the 
energy of such a state crosses the chemical potential fi 
(which is the same in all reservoirs), the QD gains or loses 
one electron, so that the total pumped charge flowing into 
it (per period) is quantized. However, the distribution of 
the pumped charge between different reservoirs is pro- 
portional to the corresponding coupling strengths (tun- 
nelling rates). Therefore the pumping current between 
any two leads can be obtained by summing up individual 
resonance contributions, with appropriate signs. 

A similar 'shuttling mechanism' for pumping has been 
used widely to interpret experiments^ in the Coulomb 
blockade regime, when the energetics on the QD is dom- 
inated by the electron-electron interactions^. In that 
approach, electrons are transferred from a lead to the 
dot and then from the dot to another lead, whenever 
such transfers are favored energetically. In contrast, 
Refs. lllUl5lllal24ll25l presented explicit quantum mechan- 
ical calculations for pumping of non-interacting electrons, 
calculated the total charge pumped during a full cycle 
and emphasized the role played by quantum interfer- 
ence in such processes. In some sense, the present pa- 
per bridges between these points of view: in the limit 
of weak coupling between the QD and the leads, we 
do end up with a loading/unloading picture, even for 
non-interacting electrons^. However, the details of the 
charge exchanges during a pumping cycle are found to 
be more complicated than in the 'shuttling' picture: at 
a given resonance point, charge can usually be shared by 
several leads. Apart from this, the conditions for the ap- 
plicability of our loading/unloading picture are similar to 
those of a single electron transistor—, in the sense that the 
role of quantum interference is restricted to the definition 
of independent single-particle resonances. In view of this, 
there is room to conjecture that some of our results may 
also apply in the presence of electron interactions. 

The paper is organized as follows. In Sec. [H] we review 



the physical assumptions of the model and the formulae 
used for the calculation of the adiabatic current. We then 
use these formulae to derive the current for a single reso- 
nant state, by approximating the Green function on the 
QD by a Breit-Wigner-type formula. In Sec. lIIII we obtain 
our main result — the adiabatically pumped charge for 
a sequence of well-defined distinct resonances — and dis- 
cuss possible applications and experimental verification. 
To demonstrate this general picture, Sec. II VI presents the 
analysis of the pumped charge for a simple model 15 of a 
"turnstile" pumping device. A short summary concludes 
the paper, in Sec. 



II. ADIABATIC CURRENT 

We consider a spatially confined nanostructure (the 
QD) connected by ideal leads to the electronic reservoirs 
with a common chemical potential /1 and temperature 
T. The total Hamiltonian for non-interacting spinless 
electrons is 

n = n d + J2{K + L a + Li) , (i) 

a 

H d = Y / h mn(t) dld n , (dot) (2) 

mn 

n l a = J2 E *k cl k c ak , (leads) (3) 

fc 

L a = X a (t) ^ Jakn c ak d n . (hopping) (4) 



Here H d is the Hamiltonian of an AT-state isolated QD 
(n, m = 1, . . . , N), the index a — 1, . . . , L enumerates the 
one-dimensional leads connected to the QD, c^ fc creates 
a standing wave \wk a ) with wavenumber fc and energy 
E a k in the channel a, the operator L a describes hopping 
from the QD into the channel a, and the A a 's are real di- 
mensionless coefficients. For pumping we allow variation 
of Tl(t) via the time-dependent parameters h mn and A a . 

The instantaneous adiabatic current in the channel a, 
directed from a remote reservoir towards the QD, has 
been expressed in Ref. as 



I a (t) 



e 

2^ 



dEf(E)l a , l a 



1 



(xk a \H\xka) , (5) 



where f(E) = 1/[1 + e (E-^)/k B T^ is the Fermi-Dirac dis- 
tribution and |xfca) is the instantaneous scattering state 
normalized to a unit flux, (xka\Xk'a) — i^/vka) S(k — 
fc'), with Vka — dEka/d(hk) being the velocity in the 
channel a. 

In Appendix A we use standard scattering theory for- 
mulae to rewrite this equation in the form 



I a (E,t) = TY d \G] 1 (H d +£)G d f a + {G d + G\) f Q /2 



(G) 
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Here, the operators 

G d =(E + iQ-H d -±)~ l , (7) 
t = i-it/2, (8) 
f = f «. f a = iit (Gj, - G^)i a , (9) 

Of 

£ = ^4, 2 4 = 4(<£ + G£)L a (10) 

a 

act only on the subspace of the QD. Also, G a denotes 
the retarded Green's function of an isolated channel, 
G l a = (E + iO - H^y 1 . We have separated the self- 
energy operator £ into a sum of resonance width and 
shift operators^, r Q and £ a , which are Hermitian. 

Equation © is a generalized version of the pump- 
ing current formula derived in Ref. 0m| for a particular 
case of single mode tight-binding (TB) leads and time- 
independent couplings. — We next use the general result 
© for the special case when the pumping is dominated 
by a single resonance on the QD. Specifically, we discuss 
localization of the accumulated charge and show that 
variations of the charge density in the leads can be ne- 
glected for a thermally broadened resonance fSec. Ill B)l . 



A. Pumping current due to a single level 

The adiabatic current © can be calculated exactly, 
provided that one is able to compute the Green's func- 
tion © on the QD. We are interested in the regime when 
the transport is dominated by a single non-degenerate or- 
bital state, as realized, for example, in a strongly-pinched 
Q D 11 ! 27 . In this case the transitions between different 
energy levels of Ti d can be completely disregarded, and 
instead of the general Eq. J2J it is sufficient to consider 



Hi s = e(t)W ^\=e(t)dU. 



(11) 



The energy distance to the next resonant state, A, will be 
assumed to be much larger than all other energies. A can 
be set either by the non-interacting effective level spac- 
ing of an arbitrary Ti. d or by the Coulomb blockade away 
from the Kondo regime 2 ^, as discussed in Ref. [n[ In the 
former case, for example, our approximation can be ap- 
plied to the much studied double-delta barrier pumping 
mo{ j e j9 i i4 i 29 - n k e resonant tunnelling regime*^. The 
energy level e(t) may be different at different stages of 
the pumping contour, as discussed later on in Sec. IIIII 

The Green's function corresponding to TC d cs now as- 
sumes the Breit-Wigner-like formSi 



Gd 



E-e - (V>|£|^) 
Substitution of Eq. ifT^) into Eq. © gives 



T 



T Q Eq — T a (Eq — E) 

(E-E )2 + (r/2r 



(12) 



(13) 



where E (E,t) =e+(ip\S\tp) and T a (E,t) = (V>|f a |V>). 
Since the partial 'width' T a is of order }? a \J a k\ 2 it rep- 
resents a measure for the coupling of the QD with the 
channel a. The exact adiabatic current for a single level 
given by Eq. l|13JI will be the starting point for our anal- 
ysis of the pumped charge in Sec. lIIII Breit-Wigner-type 
expressions for the current pumped by a single orbital 
level have been derived previously in the weak pumping 
limit^i, and in the presence of interactions and Zeeman 
splitting^!. However, they were not used to discuss the 
details of the pumped charge quantization. 

In the remainder of this section we discuss the physical 
interpretation of Eq. (|13|) . The total current, / = ^ Q I a , 
represents changes in the total charge accumulated both 
on the dot and in the leads. For small dot-lead couplings, 
one would expect that the charge on the QD itself is a 
well-defined quantity and a simple picture of single elec- 
trons tunnelling between the leads and the QD should ap- 
ply. In order to clarify the relation between our quantum 
calculation and this 'classical shuttling picture', we first 
discuss the localization of the charge and then ( Sec. Ill B"|) 
use a standard master equation approach to re-derive the 
current formula in the limit of T <C kT. 

Explicitly, Eq. Ijl3(l implies that the total current in 
our model is a full time derivative, I = dQ F (Eo, T)/dt of 
some time-dependent charge Q F (t), where 



'(*) 



./£/'(/:; ^ + larctan 2(jE r jEo) 



(14) 



(We have chosen the integration constant such that Q F /e 
is bounded between and 1.) The charge Q F can be 



mterpretec 



as the additional charge induced in the 



system by an extra impurity state \ip), in the same way 
as in the Friedel sum rule for the Anderson impurity 
model 31 i 32 . 

This delocalized charge, Q F , is to be compared 
with the local equilibrium occupation inside the QD, 
which is given by Q occ /e — Tr[p \tp) (ip\], where p = 
h^ 1 J dE f(E) J2 a \Xka) (Xka\ is the equilib rium d ensity 
matrix corresponding to 7i(t). 43 Using Eqs. (jA.lfl . (|A.6(1 . 
and l|12f> one can show that 

Q ^ = ±J dEm —± T — 2 . (15) 

If E Q and T were independent of E, then integration by 
parts would yield the equality Q occ — Q F . In general, E 
and r do depend on E, and hence Q occ ^ Q F . 



B. Alternative derivation in the limit T kT 

In the limit T <C kT, the interference effects in the 
leads are expected to be irrelevan t 23 ' 27 and the resonant 
energy level E (t) is characterized by its non-stationary 
occupation probability P(£), which changes due to tun- 
nelling. Following Ref. Li'l we assume energy conserva- 
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tion in the tunnel process and write down the total tun- 
nelling current in the direction from the lead a into the 
QD as 

I a = e(T a /h) {[1 - P]f(E ) -P[l- f(E )]} . (16) 

The tunnelling rates T a /h should be calculated at E = 
Eq. 

In equilibrium, the occupation probability P cq is 
given by the Gibbs distribution in the grand canonical 
ensemble^. For a single energy level it is the same as 
the Fermi-Dirac distribution, P cq = f(E ). Substitution 
of P = P cq into Eq. l|ltjfl gives a vanishing net current, 
as anticipated from the detailed balance. 

The kinetic equation for the non-equilibrium distribu- 
tion it; 27 - 33 



P = J2l a /e=(T/h)[f(E )-P}. 



(17) 



For slowly varying T, Eo, the time-dependent probability 
P will relax quickly towards P cq . To get an adiabatic 
solution, we substitute P = P cq + P into Eq. 117|) and 
set dP/dt = 0. The result is 



P 



-hf(E )E /r. 



(18) 



The current l|16|) corresponding to the adiabatic solution 
(HBJ) is 



I a = e(T a /T)f(E )E = e(T a /T)P ( 



eq ■ 



(19) 



(Non-adiabatic corrections can be also included, using 
the exact solution of Eq. ((17^ 1. 

Now we can compare the current derived for kT <C T, 
Eq. I|19l) . to the general expression, Eq. (|13l) . which is 
valid for arbitrary kT and T. By considering T/kT as a 
small parameter, one can perform the energy integration 
in Eq. © with I a given by Eq. (|13fl to obtain 



I a = e(T a /T)f'(E ) E + E 0(T/kT) 



f a 0{Y/kT) 

(20) 



which is consistent with Eq. (|19f) . The smallness of the 
term with r Q comes from the fact that it is multiplied 
in Eq. (|13|1 by an odd function of (Eq — E). Similarly, 
the integrals in Eqs. (|14fl and l|15|) give the leading order 
equalities Q occ = Q F = eP cq . 

We conclude that in the case of a thermally broadened 
resonance, the simple tunnelling picture which considers 
a well-defined charge on the QD is consistent with the 
general scattering approach leading to Eq. (|13fl . 



III. RESONANCE APPROXIMATION 

The Breit-Wigner form (|13|l of the pumping current 
demonstrates a well-established fac*2iiiii&, that pump- 
ing is greatly enhanced near a resonance. The resonance 



condition is \Eq — /i| < D, where D = max(r, kT) is 
the energetic width of the resonance. One option, con- 
sidered in Ref. I^jj, is to design the pumping contour in 
such a way that the system stays entirely at resonant 
transmission. In this case, the Breit-Wigner approxima- 
tion does not lead to any pumped charge quantization 20 . 
Here we focus on a more generic case, when the resonance 
condition is satisfied only during a small fraction of the 
pumping cycle, as the sys tem go es through a resonance 
point. As shown in Refs. IIIII6IITIT9I this situation al- 
lows for pumped charge quantization. Specifically, we 
assume that the system remains near a resonance point 
only during a small fraction of the pumping cycle. This 
requires relatively narrow resonances, i. e. small widths 
D and therefore also small T. 

Consider a resonance time tn on the pumping contour, 
identified by the resonance condition Eo(fi, tp) — /z. This 
identifies a 'resonance point' on the contour. Assume also 
that the system 'crosses' this resonance point completely 
between times t\ and ti, such that 

1. r a , Eq are energy independent around the Fermi 
surface (for \E — [i\ < kT); 

2. at the 'boundary' times, the system is far from the 
resonance, D < \Eq(jj,, £1,2) - Ml < A; 

3. while at resonance, the couplings change negligibly, 

|r Q |«|p |. 

Under these conditions, we can integrate Eq. I|13|) and get 
the charge transferred from the reservoir a in a simple 
form: 



dt I a = -e-p sgn E , (at E = fi). 



(21) 



For this particular resonance point, other parts of the 
pumping contour contribute negligibly to this charge. 
Equation (|21|l is our main result for the pumped charge 
due to a well-defined resonance point. We will refer to 
this result as 'the resonance approximation'. In this ap- 
proximation, each reservoir contributes on average a frac- 
tion of the electronic charge, which is proportional to the 
corresponding fractional decay width or coupling r Q /T. 
The total change in the charge accumulated in the system 
due to this particular resonance is thus 



±e. 



(22) 



This result can be easily generalized for several inde- 
pendent resonance points. If the pumping contour can 
be separated into several parts, each containing a single 
well-defined resonance point, and if the pumping currents 
on the rest of the contour remain negligible, then the to- 
tal charge Q^, pumped through the channel a, is given 
by a sum over the resonances: = X^res ^Qa S - For a 
periodic Ti.(t), the pumping contour is closed, and charge 
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conservation J2 a = is ensured by Eq. ({2*2*)) and the 
fact that the number of loading (Eq < 0) and unloading 
(E > 0) resonance points is the same. 

A. Pumped charge quantization 

Equation (|22[1 can be interpreted as the load- 
ing/unloading of exactly one electron into/out of the QD, 
depending on the sign of Eq at the Fermi level. Further- 
more, Eq. (|21|l implies that AQ rcs is dominated by the 
current from a single channel a, provided that Y a 3> T Q ' 
for a a'. If the same applies to all the resonances, 
then we end up with a "classical" picture, in which the 
pumping cycle contains a sequence of individual discrete 
events, of exchanging electrons one by one between a 
reservoir and the QD. After a full cycle, the charge on 
the QD will remain unchanged, and an integer number 
of electrons will have crossed the QD between any pair 
of reservoirs. This gives a detailed explanation of the 
pumped charge quantization within this approximation. 

Using the same conditions as used to derive Eq. (|21|l . 
one can show that both AQ F = Q F (t2) — Q F (t\) and 
AQ occ are equal to AQ ros . This means that every time 
the system crosses a resonance point, the charge associ- 
ated with the resonant state changes by ~ ±e. There- 
fore we stress that if one is interested in the total charge 
pumped by a single resonance (and not, for example, in 
the lineshape of the current, Eq. then the simple 

picture of loading/unloading of a single electron, as re- 
flected in Eq. I|21|) , is applicable — regardless of the ratio 
T/kT. 

We also note that for such an ideal quantization (Q^ — > 
e x integer), that is independent of the contour details, 
one would need to consider the limit r a i=s T — > for 
each resonance; the resonance approximation becomes 
exact, with results which are independent of the details 
of the contour, when r — > 0, and the charge goes only via 
channel a when r Q /r — > 1. As explained in the next sub- 
section, this implies a vanishing transmission throughout 
the whole pumping cycle, in accordance with the conclu- 
sions of Refs. Illl35l 



B. Relation to conductance 

The criteria for the validity of the resonance approx- 
imation, listed in the previous subsection, can be quan- 
titatively checked in experiments (or in numerical calcu- 
lations) by monitoring the conductance between differ- 
ent leads as a function of parameters along the pumping 
contou»iiii&. A definitive signature of the relevant trans- 
port regime (for having a significant non-zero pumped 
charge) would be the presence of an even number of well- 
separated peaks in the conductance time trace: each res- 
onance (M in Fig. ["J is associated with two peaks in 
the instantaneous transmission, encountered at the two 
resonance points (B and D) where the pumping contour 



crosses the resonance line on each side of the resonance, 
as schematically shown in Fig. ["] Note that this mea- 
surement is independent of time: one simply measures 
the conductance at different points on the pumping con- 
tour. 

The contribution of each particular conductance peak 
to the pumped charge can be calculated along the fol- 
lowing lines. Application of the general expression of the 
transmission probability^ from channel a' to channel a, 
T aa , = - J dE f'(E)Tr[Glt a G d T a ,], to our resonance 
model (as defined in Eq. gives the standard Breit- 

Wigner^S result (see e.g. Ref. l27j) : 



'EqlQl' 



dEf'(E) 



(E - EqY + (T/2f 



(23) 



Let us consider for simplicity an example of L single- 
mode leads. By using the multi-terminal Landaucr 
conductance formula^ for spinless electrons, Gaa' = 
(e 2 /h)T aa ', in Eq. (|23fl . we recover well-established^ re- 
sults for the peak conductance of a strongly pinched QD, 
that are related to Eq. (|21|l in an extremely simple way: 



■>pcak 



e 2 4F F / 
~h YD 



AY 



AQ™ S AQ 



res 
a' i 



where 



' (8/ir)kT, fcT>r. 



(24) 



(25) 



Measurements of the peak conductance at a particular 
resonance point for fixed temperature and all possible 
combinations of source and drain leads would give, in 
principle, (L 2 — L)/2 experimental values to be used in 
Eqs. C2J). Together with Eq. ("""""J, this gives (L 2 -L)/2 + 
1 equations for the L + 1 unknowns AQ^ CS and Y/D. 
Measurement of the temperature dependence of G^Li (T) 
would yield D(T), and thus determine Y. We see that 
even for L = 2 it is possible to predict the adiabatically 
pumped charge from the conductance measurements, and 
for L > 2 different cross-checks become feasible. 

Additional input of a few bits of information is neces- 
sary to make the solution of Eqs. I|24() and (|22() unique. 
For a specific resonance 'res', all the charges AQ™ S (for 
all a) have the same sign, determined by the type of the 
resonance: "+" for loading and "— " for unloading, see 
Eq. H21[l. An additional sign uncertainty arises in the 
case of two terminals (a — l,r): the respective equation 
for the pumped charge, AQ t {e - AQ r ) = Qf^(hD/AY), 
is symmetric under inversion, I r. The resolution of 
these uncertainties depends on the particular experimen- 
tal situation, and should be easy in simple cases. We 
illustrate this point in Sec. II VI below, when we discuss a 
two-terminal example. 



C. Adiabaticity condition 

One condition for the validity of the adiabatic picture 
requires that an electron should have enough time to tun- 
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nel under the barriers while the system is at resonance. 
Thus, the inverse tunnelling rate h/Y should be much 
smaller than the duration of the resonance, r r = D/\Eq\, 
yielding the adiabaticity condition, 

h\E \^YD. (26) 

This condition implies that both the amplitude and the 
frequency of the pumping potential must be sufficiently 
small for an adiabatic pump2&. The resonance duration 
r r can be extracted from measurements of the conduc- 
tance as follows: measuring the variation of the con- 
ductance through the resonance, using a very low fre- 
quency ojq, would yield the resonance width r r o for that 
frequency. The value of r r relevant for the pumping ex- 
periments can then be found by rescaling, r r — t t qujq/uj. 

At zero temperature, D = Y and the condition l|26() 
can be compared to the adiabaticity criterion for co- 
herent pumping formulated recently by Moskalets and 
Biittiker 29 . They consider the number of side-bands n maK 
required to describe adequately the Fourier transform of 
the instantaneous scattering matrix. In our case the res- 
onant peak of transmission in the time domain has the 
width t>, and the number of relevant Fourier harmonics 
"max is at least (wT r ) _1 , where lu is the cyclic freq uency 
of the pump. The adiabaticity criterion of Ref. |29| states 
that the scattering matrix should vary little with energy 
over the range E ± /zwn max . Since our characteristic en- 
ergy scale for the scattering matrix is Y, the condition of 
Ref. l29l takes the form Y 3> Sum max = Kt~ , equivalent 
to Eq. ijSBT) . 

For a thermally broadened resonance, the adiabaticity 
condition follows also from the requirement that the first 
order correction P to the adiabatic solution of the kinetic 
equation (JT7|l remains small: P < 1 \f'(E )E Q \ < 

Y/h^> Eq. pel). 



the Fermi energy \x. Therefore, the following algorithm 
can be formulated: 

1. Diagonalize Ti d (t) (analytically or numerically), to 
get the spectrum {e m (i), \ip m (t))}. 

2. Calculate the time-dependent decay widths 
r?(t) (ip m (t)\Y a (E = n)\i> m (t)) and 
shifted energy levels e' m (t) — e m {t) + 
(<ip m (t)\E(E = fM)\i> m (t)). 

3. For every m, find all such times i mj when the res- 
onance condition e' m (t m j) = fi is satisfied. 

4. At each resonance time t — t m j, compute the cor- 
responding partial charge <?™ J = er™/^ a , T™,. 

5. Calculate the total pumped charge as 

Qa=-EC Js g n 4(< mJ ), (27) 

or set = if no resonances were found in step 3. 

The application of this algorithm is justified under the 
conditions listed in the beginning of this section. The 
most important condition is the consistency of the per- 
turbation expansion, Y™(t m j) <g; A(t mt j), where A(t) is 
the level spacing oiTL d {t) at the Fermi surface. 

The algorithm will fail for certain values of the ad- 
justable (not pumping) parameters of the model, for 
which the number of resonance points found in step 3 
changes. This change corresponds to the appearance (or 
annihilation) of a pair of loading/unloading resonances. 
Such a crossover is usually manifested by a sharp change 
(a step)i§i2i in the total pumped charge, as function of 
the model parameters. 



D. Application to complicated pumping potentials 

In the resonance approximation, the pumped charge 
is expressed in terms of the resonance points, where the 
pumping contour crosses the resonance lines, and do not 
require the full information on the contour in the param- 
eter space. We now discuss the conditions under which 
Eq. i|21|) can be used to obtain efficient approximate esti- 
mates of the pumped charge for a model Hamiltonian 
Tt d , which is complicated enough to render an exact 
integration^ of Eq. @ impractical. Even when the valid- 
ity of the resonance approximation is marginal, such an 
approximate estimate could provide a handy t ool for ex- 
ploring complicated pumping models (e.g. Ref . GHHHi 
and identifying the relevant physical parameters. For 
simplicity, we restrict this discussion to zero tempera- 
ture. 

To leading order in the coupling strengths A Q , the pa- 
rameters of the resonant level in Eq. are given by the 
eigenstate of the decoupled Ti. d (t) which is the closest to 



IV. EXAMPLE: TURNSTILE MODEL 

We illustrate the resonance approximation by a simple 
example of a single energy level with adiabatically vary- 
ing couplings to the left and right reservoirs (single level 
turnstile modeli^). Applications to more complicated 
models, such as pumping by surface acoustic waves^, 
will be reported elsewhere. 



A. The turnstile pumping model 

The single level turnstile model, discussed in Ref. 0, 
can be described as a special case of the general Hamil- 
tonian with N = 1 site (and a single energy /in = e) 
on the QD and with L = 2 leads, denoted by a = l,r. 
It is now convenient to use a slightly different notation: 
Consider an infinite chain of tight binding (TB) sites, 
enumerated by n — 0, ±1,.... The site n = 0, which 
represents the QD, has a time-independent energy e and 
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defines Ttf cs = ed^d, with eigenstate \ip). The sites with 
n > (n < 0) form the right (left) single-mode TB lead: 

(<4c„±i + c, T i±1 c„) , (28) 



H„ — 



E 

n=±l 



J I 



where the upper sign refers to a = r. The coupling 
operators are L Q = y/X\JJ)J c'^d, with the two time- 
dependent pumping parameters X a — A„ . 

The Hamiltonian of the leads l|28|) is characterized by 
the dispersion relation Ek = —2 J cos ka and the retarded 
Green's function 



[G l a ] 



Aka\n— ml 



u ika\m-\-n\ 



(29) 



i2Jsinfca 

where a is the nearest-neighbor distance. The self-energy 
operator [Eq. ©] is t = -{X t + X r ) J e lka dU. 

We consider the zero temperature limit and 
parametrize the on-site energy as e = (— 2 + 8) J cos Ka, 
where the dimensionless parameter <5 is a measure of the 
de-tuning of the isolated level e from the Fermi energy 
/.t = — 2 J cos Ka in the leads. Near the band bottom one 
has 8 ~ (e — n)/ J. 

The resonance parameters at the Fermi surface are 

T Q = 2X a J sin Ka , 

E = {-2 + 5- Xi - X r ) J cos Ka. (30) 

Both r a and Eq depend on time via the time-dependent 
couplings X a , which span the parameter space {Xi, X r }. 
The resonance condition E = fj, defines the resonance 
line, Xi + X r = 8. For an explicit calculation, we next 
choose the pumping contour to be a square with corners 
at points A{Xi;Xi) and C{X2] X2), as used in Ref. ITEI 
(This is shown in Fig. [2k-, which forms an explicit example 
ofFig.UJ. 

The necessary conditions of Sec. Mil for having distinct 
resonances, are satisfied only at the bottom of the TB 
band (sin Ka <C cos Ka) . As we gradually increase 8 from 
zero, the resonance line in Fig.|5Ji moves in the direction 
indicated by the small arrow. The resonance line crosses 
the contour only if 2Xi = 5i < 6 < 83 = 2X 2 . Therefore, 
within the resonance approximation we have 

Q R /e = 0, if 5 < 81 or 5 > S 3 . (31) 

For the direction of the contour shown by the arrows 
in Fig. [21 the resonance point B corresponds to loading 
of the dot mostly from the left (r^ > T r ). Its comple- 
mentary resonance point D is associated with unloading 
mostly to the right (r^ < T r ). This interpretation is 
illustrated in Fig. |2Jd. 

At the lower left part of the contour, 5 < 82 — {5\ + 
83)/2 = X\ +X2, the resonance points are D{X\, 8 — X\) 
and B{8 — X\,X\). The partial charges pumped from 
the left [using Eqs. (2U and El)] are 



AQf = -e 



2Xi J sin Ka _ _ Xi 
2{Xi+X r ) J sin ku & 8 ' [ ' 



AQf = e- 



Xi 



(33) 



mX,. 



x 2 - 



Pumping contour 
Resonance line 



H *■ 



(a) X^ 



X 2 



E (t) 



(b) -2J 



fi=E {U) 




\v=Eo(t-)\ 



FIG. 2: (color online) (a) The pumping contour A-B-C-D- 
A and the resonance line B-D for the single level turnstile 
model^. (b) Interpretation of the pumping cycle on an en- 
ergy diagram. A: The effective energy level Eo is above the 
chemical potential fi, the dot is empty. B: Loading process 
with preference to the left-coming electrons. C: The level 
Eo is below /i, the dot is occupied. D: Unloading process 
with preference to the right-going electrons. The asymmetry 
between B and D creates the non-vanishing total pumped 
charge. The arrows indicate schematically the direction and 
the relative magnitude of the current pulses caused by each 
resonance. 



where we have used sgn Eq = sgn ^ {—Xi — X r ) = +1 for 
point D. The net pumped charge is thus 

Q R /e = Qf/e = (AQf + AQf )/e - -Q?/e 

= l-{8 1 /5), if 8!<S< 8 2 . (34) 

A similar analysis for crossing at D{8 — X 2 ,X2) and 
B{X 2 , 8 - X 2 ) (when S 2 < 8 < S 3 ) yields 



) K /e = {8 3 /8) - 1, if <5 2 < 5 < ^ 3 



(35) 



Our resonance approximation results for Q R /e are 
shown for some typical parameters (together with the 
exact results, see below) in Fi g. EH These results agree 
qualitatively with those of Refs Illl andlT^: Q R /e reaches 
its maximum value {X2— X±)/{Xi+X2) at 8 = 82, where 
the resonance points B and D are farthest away from 
the resonance point M, which occurs at X\ = X2 = 8/2. 
Note that Q R /e approaches the quantized value 1 when 
X 2 /Xi — > 00, i.e. when the transmission at the resonance 
points (related to AX 1 X 2 /{X 1 + X 2 ) 2 , via Eq. (J2SJ) van- 
ishes. This is consistent with Ref. 11, which required that 
"a large part of the resonance line" be surrounded by the 
pumping contour. 



Q 



1 - 
X2-X1 

X 2 +X! 








§1 §2 83 



FIG. 3: (color online) Pumped charge (in units of e) as a 
function of 8 for X\ — 1/50, X2 = 1/5 and na — it/20, calcu- 
lated within the resonance approximation (Q R , blue contin- 
uous line) and exactly (Q, dashed line). Thick bars on the 
8 axis mark the resonance widths izT/J around the special 
points S\ t 2,3, where deviations from the the exact result are 
anticipated. Inset: Absolute error of the resonance approx- 
imation, (Q R — Q)/e, for the same values of 5. The thick 
dotted line corresponds to Q R calculated from the transmis- 
sion maxima, see text for details. 



8 




4 1 4 1 • v 

A C A g 
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FIG. 4: Time traces of the transmission coefficient % r along 
the pumping contour for six values of 5, increasing with con- 
stant intervals from top to bottom. Two complementary 
resonances B and D (marked with arrows) are observed for 
5\ < 5 < 5z when the pumping contour crosses the resonance 
line (cf. Fig. Ei). Thick bars on the 5-axis mark the regions 
where the loading/unloading pumping mechanism fails. 



B. Comparison with exact results 

The formula (|13|) for the resonance current is exact 
in our case. Substitution of Eq. I|3l)[l into Eq. <|13|) and 

integration over the contour A-B-C-D gives the total 
pumped charge in the form 



Q = ^JdX [F(X, Xx) - F(X, X 2 



(36) 



where 



F(X, Z) = 



(8 — 2 Z) sin na cos na 



(8 - X - Z) 2 cos 2 na + (X + Z) 2 sin 2 na 



(37) 

This result was obtained in Ref. [15| using the time- 
derivatives of the scattering matrix. 

In Fig. we compare the exact Q and the approximate 
Q K . As the resonance line in Fig. [5] moves from point A 
to C, the pumped charge rises from zero to a maximum, 
close to (X2 — X\ )/{Xi +X\ ) , and then falls back towards 
zero. Except for the vicinity of the special points 8 — 
81, 82 and (S3, there is an excellent agreement between 
Eqs. and Eq. 

The most significant source for deviations of the exact 
pumped charge Q from the separated resonance result 
Q R is the term proportional to r Q in the expression of 
the pumping current (|13f) : 



e 

2^ 



dV a (E - n) 
(p E Q ) 2 + (T/2)* 



(38) 



In our example, Eq. (|5U|> yields dT a /dE — T a /E Q = 
— 2tan(«;a) when the resonance is on an edge of the con- 
tour curve where only X a varies. Thus, the integral in 



Eq. is negligibly small as long as the distance be- 

tween the resonance point and a corner of the square con- 
tour is larger than T/J. Indeed, this agrees with Fig. EI 
where the regions \8 — 8i\ < T/J are indicated by hori- 
zontal bars on the 5-axis. 

Figure 01 shows the Breit-Wigner transmission coeffi- 
cient Ti r , calculated from Eqs. (|23|) and l|3U|) as a func- 
tion of time (defined homogeneously along the pumping 
contour) for several values of 8. As 8 is increased from 
zero, a single peak develops at 8 — 81, then splits into 
two independent resonances B and D, which move along 
the pumping contour and finally merge at 8 = 83 and dis- 
appear. By comparing Fig.^to Fig.[3]one can follow the 
correlation between the presence of separate well-defined 
transmission peaks and the validity of the resonance ap- 
proximation for the pumped charge. 



C. Relation to transmission 

The quantitative relation between the pumped charge 
and the transmission (conductance) has been discussed 
in Sec. 1111 Bl In order to illustrate this discussion, let 
us assume that the transmission traces (Fig. are the 
only available data for our two terminal system. One 
observes two resonances in the range 8\ < 8 < 83 - both 
giving the same value of the peak transmission T max . 
One of the resonances represents loading, contributing 
AQ^ CS > 0, while the other one necessarily represents un- 
loading (with AQ r ° s < 0). If we make a mistake at this 
stage and take the wrong sign in Eq. I|22|) . it will only 
change the assumed pumping direction, Q R — > — Q R . 
Let us treat the first resonance as loading and calcu- 
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late the partial charge pumped from the left reservoir, 
AQ[ es > 0. Solution of Eqs. g2J| and O gives two 
roots, AQ rcs = e(l ± y/1 — 7^ nax )/2, and one must de- 
cide which of the two corresponds to AQ[ es . The same 
dilemma holds for the second resonance. Considering all 
four options yields three possible answers Q R = ±Q' and 
0, where Q 1 = e\J\ — T max . The correct result (Q K = Q') 
may be chosen as the one which gives the best fit to the 
data of the pumping calculation/experiment. Once the 
uncertain signs have been chosen correctly, there is no 
need to repeat this "trial-and-error" procedure, since the 
contour changes continuously. Of course, if some fea- 
tures of the pumping contour design are known (such as 
which coupling is dominant in different regions), the sign 
uncertainties are much easier to resolve. 

The result of the above calculation, \Q — Q R \/e, is 
shown in the inset of Fig. by a thick dotted line. One 
can see that both ways of calculating Q R (from the an- 
alytic expressions (|3 1134135ft and from using the peak 
transmission) give similar small deviations from the exact 
value Q of the pumped charge. 

We now leave our specific example, and consider Eq. 
(|38[l for a general resonance. As seen in the example, 
the integral in Eq. I|38[) becomes non-zero whenever 
7a = r Q /i?o is not a time- independent constant during 
the whole resonance. For non-constant j a , the largest 
deviation \Q — Q R \ arises when 7 a changes sign exactly 
at the resonance point Eq = fx; one then finds that 
\(Q — Q R )/Q R | < xmax|7 Q |, where x is a number of 
order unity, which depends on the details of the contour. 
These considerations justify Condition No. 3 in the be- 
ginning of Sec. IIIII 



V. CONCLUSIONS 

We have considered a general model of adiabatic quan- 
tum pumping of spinless non-interacting electrons, in the 
coherent resonant tunnelling regime. In the limit of dis- 
tinct transmission resonances along the pumping con- 
tour, the pumped charge is given by a sum of individual 
contributions due to each resonance. During each res- 
onance one electron either enters or leaves the system, 
with the probability distribution between different reser- 
voirs given by the corresponding tunnelling rates T a /h. 

We have clarified the role of quantum coherence in the 
resonance-assisted pumped charge quantization by show- 
ing that quantization arises due to population of discrete 
resonant states with preference to a single reservoir in 
each resonance. A quantitative and experimentally veri- 
fiable relation between the pumped charge and the peak 
conductance has been proposed. The resonance approxi- 
mation also provides a simple calculational algorithm for 
analyzing complex pumping potentials. 

Our results remain valid if (1) the spacing A between 
different resonant levels is much larger than r a , kT; (2) 
the relative magnitude of the couplings T a to different 
reservoirs does not change much during a resonance; (3) 



the adibaticity condition HEq <C max(r, kT) is not vio- 
lated. 

Systematic extension of the resonance approximation 
to the situations when electron-electron interactions can 
not be neglected is a question of further study. 
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APPENDIX: FORMULA FOR THE PUMPED 
CURRENT 

In this Appendix, we use standard scattering theory 
relationa 26 ! 37 to derive Eq. JBJ from Eq. The scat- 
tering states \xka) can be obtained from the Lippman- 
Schwinger equation 



\x ka ) = {l + GLl)\w ka ) 



(A.l) 



where G = (E+iO—Tt) -1 is the retarded Green's function 
taken at energy E — E ka . The time t enters Eq. i A. 1|> 
as a parameter. 

Defining projection operators P d and P a onto the QD 
and onto lead a, one has 

U d = pd U dpd^ Gd = pd Gd p d ^ p aU p ai = Saa ,n l a , 

L a = L a P d 
and therefore 



P a HP d , P d \w ka ) = 0, (A.2) 



M d a = (xk a \n d \ X k a ) = {w ka \L a G d n d G d Ll\w ka ). 

(A.3) 

To derive Eq. (7), we start from (E - H)G = I, mul- 
tiply from the right by P d and from the left by P a and - 
using the identity P d + -P a = I - obtain the relation 
P a GP d = G l a L a Gd- A similar multiplication from the 
left by P d then yields Eq. Q, with 



S — ^ L\G l a L a , 



(A.4) 



which is equivalent to Eq. (|SJ. 

Similarly, the time dependance of the coupling 
strengths X a (t) contributes to the current I a via the 
matrix element M l a = J2 a >(xka\L a >\xka) + h.c. Using 
the trivial relations L a — (X a / 'X a )L a and P a > \Wk a ) 



■ 



a straightforward calculation gives: 



ML 



(Wkc 



G d £G d Lt 



(G d + G d ) 



(A.5) 



The normalization to the unit flux, 
(wka\wk' a ) = (27r/ufc a ) 5(k — k') implies that 
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P a = J(dk/2Tr)vka \wka) {wka\- Using also the standard 
relation i(G l a — G^J) = \wk a ){wka\/K we find the 
relation 

ht a = L a \w ka ) (w ka | L a . (A. 6) 



Introducing the trace over the QD's subspace, using (|A. 6|) 
in Eqs. I|A.3I) and (|A.5|I . and substituting the results into 
J a = (M* + M l a )/h, we finally end up with Eq. ® . 
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Distinct charge steps during the pumping cycle have been 
presented in Ref.l24l but with no detailed quantitative dis- 
cussion of their relation to resonances. 
Eq. JSJ can be used to demonstrate the equivalence be- 
tween the formalism used in Refs. Il5lll6u24ll25l and the adi- 
abatic limit of the Keldysh nonequilibrium Green's func- 
tion techniques—*—. Without loss of generality one can 
set X a — 0, and use the cyclic property of the trace in 
Eq. JSJ. Extending the trace to the whole system and mul- 
tiplying the result by two to restore spin-degeneracy gives 
the charge injected into the system during an infinitesi- 
mal time St as 8Q a = -(e/vr) / dEf (£)Tr[G f f a GH5i\. 
This is identical to Eq. (4) of Ref. 123 . which in turn was 
proved—*— equivalent to the Brouwer formula**— 
The use of the equilibrium p is justified because adiabatic 
pumping is invariant to the rescaling of time and can al- 
ways be performed quasi-statically. 



